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ABSTRACT 

We present three bright z + dropout candidates selected from deep Near-Infrared (NIR) imaging of 
the COSMOS 2 square degree field. All three objects match the 0.8 — 8^m colors of other published 
z > 7 candidates but are three magnitudes brighter, facilitating further study. Deep spectroscopy of 
two of the candidates covering 0.64— 1.02/um with Keck-DEIMOS and all three covering 0.94-1.10^m 
and 1.52-1.80/xm with Keck-NIRSPEC detects weak spectral features tentatively identified as Ly-a at 
z = 6.95 and z — 7.69 in two of the objects. The third object is placed at z <~ 1.6 based on a 24/im and 
weak optical detection. A comparison with the spectral energy distributions of known z < 7 galaxies, 
including objects with strong spectral lines, large extinction, and large systematic uncertainties in 
the photometry yields no objects with similar colors. However, the A > lytmi properties of all three 
objects can be matched to optically detected sources with photometric redshifts at z <~ 1.8, so the 
non-detection in the i + and z + bands are the primary factors which favors a z > 7 solution. If any 
of these objects are at z <~ 7 the bright end of the luminosity function is significantly higher at z > 7 
than suggested by previous studies, but consistent within the statistical uncertainty and the dark 
matter halo distribution. If these objects are at low redshift, the Lyman-Break selection must be 
contaminated by a previously unknown population of low redshift objects with very strong breaks in 
their broad band spectral energy distributions and blue NIR colors. The implications of this result 
on luminosity function evolution at high redshift is discussed. We show that the primary limitation 
of z > 7 galaxy searches with broad filters is the depth of the available optical data. 
Subject headings: Galaxies: Formation, Galaxies: High-Redshift, Galaxies:Evolution 
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1. INTRODUCTION 

The z > 7 universe holds the keys to many open ques- 
tions about galaxy formation and evolution. Cold Dark 
Matter (CDM) models predict that gal axies assemble 
through merging <>f dark ma tter haloes (ppringel et a I 



20051 iCole et al. |2000| iKauffmann et al. 1993} and the 



distribution of initial halo masses is strongly constrained 
by the Cosmic Microwave B ackground (CMB) measure- 
ments (Bennett et al. 2003). However, populating the 
dark matter haios with galaxies is done in a semi-analytic 
manner with many free parameters that are tuned to ob- 
servations. At z ~ 7 the universe was less than a billion 
years old, leaving little time for galaxies to evolve, so the 
luminosities and masses of galaxies at z > 7 are a direct 
probe of the assumptions underlying galaxy formation 
scenarios in structure formation models. 

Advances in red-sensitive Charge-Coupled-Devices 
(CCD s) have allowed studies of the universe out to z ~ 7 
(e.g. |Ouchi et al.|[2009{ |Ota et al.|[2008j |Iye et al.|[2008| ; 
however, accessing higher redshiits is difficult due to 
the lack of wide field near-infrared cameras and high 
sky backgrounds beyond 0.8//m in the case of ground 
based telescopes. Ultra-deep studies from space with 
Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS) and Advanced Camera for Surveys (ACS) 
on the Hubble Space Telescope (HST) h ave constrained 
the faint-end of the luminosity function (Bouwens et al 



2008| Oesch et al.||2010 ) with deep ground based studies 
and NlCMUS parallel surveys placing useful limits on its 



bright end ()Stanway et al.||2008b} |Mannucci et al.||2007 



Henry et al.||2009" ), and tensing surv eys also providing 
useful constraints on fainter galaxies flKneib et al.||2004 
Richard et aI1p008l |Bradley et aT]|2O0H| I Bouwens et, al | 
2009p albeit with significant uncertainty due to the small 
areas covered. Unfortunately, the selection function is 
complicated due to contamination from faint galaxies at 
1.4 < z < 2.5 where the galaxy population is also poorly 
understood and the fact that no objects have b een spec- 
troscopically confirmed ( Stanway et al.|[2008b ). 

In this paper we use deep IN IK imaging of the 2 square 
degree Cosmic Evolution Survey or COSMOS field with 
the Wide-Field Infrared Camera (WFCAM) instrument 
on United Kingdom Infrared Telescope (UKIRT) and 
the Wide-Field Infrared Camera (WIRCAM) instrument 
on the Canada- France-Hawaii- Telescope (CFHT) along 
with deep Spitzer Infrared Array Camera (IRAC) and 
and Multi-Band Imaging Photometer (MIPS) imaging 
to color select candidate z-dropout galaxies. We find 
three candidates brighter than Jab < 23.5, which have 
0.9-8/im SEDs consistent with fainter z > 7 c andidates 
discovered in the GOODS and UDF fields (|Bouwens| 

2008b[ |Oesch et al.||2010j 



et al. 2008 Stanway et al.| 



and have mid- infrare d (MIR) colors consistent with the 
selection proposed in 
present deep spectra o 
DEIMOS in the 0.64 



Stanway et al. 



_ ([2008bj. We then 
these objects obtained with Keck- 
1.02/im range and Keck-NIRSPEC 
in the 0.94-1.10/xm and 1.52-1.80/im range. In addition 
to these three candidates we find over 3500 objects which 
are at z ~ 1.8 and would be selected as z > 7 objects 

Associated Universities, Inc ; the 30 m telescope of the Institute 
for Radio- astronomy at Millimeter Wavelengths ( IRAM ), 
which is funded by the German Max-Planck-Society, the French 
CNRS, and the Spanish National Geographical Institute. 
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Fig. 1. — The infrared colors of our z > 7 candidates are shown 
with objects from the literature and foreground objects. In the 
top pannel the J — H and z+ — J col ors of our can didates (red 



points ) along with |Bouwens et a.1. | (|2008p (black points), |Oesch et al.| 
j |2009| l( blue points) and WFU-LlbF |Oesch et al. (2010) (magenta 
poin ts) are shown . The redder colors m |fJouwens et ah] ( |2008| l 
and |Oesch et aL] ( |2009[) are due to the use ot FllUYV, which is 
considerably wider and bluer than the J and F125W band used in 
this work and the WFC3-UDF respectively. The colors of optically 
detected sources at z < 6 with > 10ct detections in J, H, and K s 
are shown as small green points. In the bottom panel J — H and 
3.6 — 4.5 col ors of our candidates (red points) and those from|Labbe| 
|et al.||2010[ | ( blue and purp l e poin ts), |Gonzalez et al.] l |2010[ | (Green 
VointsJ, and iLabbe et al.l J2006| (black points) are shown. The 
objects presented here have colors consistent with those selected in 
deeper surveys. 



TABLE 1 
Positions 



ID 


RA(J2000) 


DEC(J2000) 


FWHM " 


Stellarity 


1 


10' 1 01 m 36.317 s 


+1:37:00.42 


2.5 


0.98 


2 


9 h 59 m 25.507 s 


+2:21:27.79 


1.5 


0.96 


3 


9 h 59 m 14.609 s 


+2:27:36.57 


2.0 


0.53 



were it not for the large difference in depth between the 
near-infrared and optical bands. Using this sample of 
interlopers, we discuss the sources of contamination for 
z > 7 studies and strategies for removing them. 

Throughout this paper we use AB magnitudes or fluxes 
in units of Jansky and a £l m = 0.3, Q v = 0.7, H Q = 70 
km s _1 Mpc -1 cosmology. 
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u* Bj g + Vj r + i + z + J H K s ch1 Ch2 




Fig. 2. — Postage stamps of the three candidate objects covering 7.5" X 7.5" in u*, Bj, g + , Vj, r + , i + , z + , J, Ks, IRAC Chi 
(3.6/xm), and IRAC Ch2 (4.5/jm) bands scaled to a common noise level and S/N stretch in all images. Notice the absence of the objects 
blueward of J band despite the significant increase in the image depth. 



2. IMAGING DATA AND OBJECT SELECTION 

The 2 square degree COSMOS field is covered by a 
combination of 30 broad, intermediate, and narrow Ul- 
traviolet (UV), Optical, and NIR filters including data 
from Galex, CFHT, Subaru, HST and Spitzer. These 
observ ations are described in detail in Zamojski et al.| 



(2007) (GALEX UV), |Capak eTal] ([20071 1 20081) (Subaru 
CTTTT, an d UKIRT optical and NIK), iKoekemoer et al. 



(12007] andlS c o yille et al.|p007| (HST optical and NIR), 



Sanders et al.| ( |20(j7 ) (Spitzer near and mid-Infrared). In 
addition there are extensive, X-r ay data obtained with 
the Ch andra X-ray observato ry (Elvis et al.| 2009 ) and 



XMM (Hasinger et al. 2007 



IRAM 30m {B crtoldi etal 



sub-mm data from the 
along with 20cm 
in 



20071 

90cm (Smolcic et al 



( jSchinnerer et al.||2007 ) and 

prep) radio images from the Very Large Array (VLA) 

Thro ughout this paper we follow the filter notation 
used in Capak et al. (20071 . In particular Subaru-Sloan 
filters are noted with the super-script + (e.g. z + ) and 
CFHT Megaprime filters are denoted with the super- 
script * (e.g. u*) to differentiate them from each other 
and the true Sloan Digital Sky Survey (SDSS) filters 
which have somewhat different bandpasses. 

This paper also makes use of new deep J band da ta ob- 
tained with the WFCAM instrument on UKIRT ( |Hirst| 
|et al.||2006[ ), along wit h H and K s data obtained with 
W1RCAM on CFHT dPuget et al.| | 2004D telescopes as 
part of the COSM OS survey ( |McCracken et al.||2009 
Capak et al. 2010[) . The J, H and K B data were taken 



over the course of three years with over 500 exposures 
at each position. This is more than sufficient to remove 
cosmic rays, moving objects, satellites, and other tran- 
sient objects that may be mistaken for high-z objects. 
Electronic cross-talk artifacts near bright objects are a 
significant problem in all three data sets. The nature 
of the cross-talk in the H and K s data allow the arti- 
facts to be modeled and removed in all but the most 
severe cases (McCracken et al. |2009|). Any remaining 



cross talk is easily identified because it consists of 64 
equally bright/dark images at 64 detector pixel intervals 
in Declination. Cross-talk artifacts are still present in 
the UKIRT WFCAM data. In WFCAM data the cross- 
talk artifacts occur at 128 detector pixel intervals in ei- 
ther Right Ascension or Declination from bright objects 
depending on position within the array and can be iden- 
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Fig. 3. — ACS F814W postage stamps of the three candidate 
objects covering 5" X 5". Note the lack of any significant detection 




Fig. 4. — MIPS 24/xm postage stamps of the three candidate 
objects covering 30" X 30". Note only object 2 is detected with a 
flux of 66 ± 12/^Jy. 



tified by their distinct positive-negative signature. 

A K s ban d selected object catalog was generated with 
SExtractor ( jBertin fc Arnouts|[l996"| using the CFHT- 
WIRCAM K s image tor detection a nd PSF matched im - 



ages for photometry as described in Capak et ah] (2007) 
The SExtracto r detec tion settings are those used m |Mc 



Cracken et al. (2009). The K s band image was chosen 
over J band tor detection because while the image depths 
are similar, the image quality is significantly better in 
the K s band images. Moreover, the expected object col- 
ors mean they should be equally bright or brighter in 
K s than J. In addition to a detection in K a band, we 
also require a 5a measurement in a 3" aperture for both 
the J and K s bands for an object to enter our analy- 
sis. This selection improves the robustness of our object 
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TABLE 2 

Photometry 



ID 


i+ 


z+ 


J 


H 


K 


m3.6 


m4.5 


m5.8 


m8.0 


1 


> 27.6 


> 25.5 


23.21 ±0.05 


23.09 ±0.12 


22.14 ± 0.02 


21.89 ±0.02 


21.65 ±0.03 






2 


26.3 ±0.5 


25.22 ± 0.5 


22.81 ± 0.1 


23.00 ±0.11 


22.11 ± 0.05 


20.96 ±0.01 


20.67 ± 0.01 


21.12 ±0.07 


20.74 ± 0.14 


3 


> 27.6 


> 25.9 


23.11 ±0.1 


> 23.5 


22.80 ± 0.14 


22.34 ± 0.03 


22.16 ±0.05 


21.9 ±0.1 


22.2 ±0.7 



Object 1 Object 2 Object 3 




Fig. 5. — The positions of the NIRSPEC (red) and DEIMOS (blue) slits are overlaid on the CFHT K s band images with North up and 
East to the left. The targeted object is marked with a red arrow. 



identification, removing a significant number of artifacts 
which are only present in either the J or K s band im- 
ages. The resulting catalog has a 5tr limiting magnitude 
in a 3" aperture of 23.7 on the AB system and is 80% 



complete a t this flux level for point sources (|McCracken 
eTat][2009). 



i'he H band data was not available when the objects 
were selected, but is used in later analysis. The pho- 
tometry was measured in the same way as the for the J 
and K s bands and has a 5a limiting magnitude in a 3" 
aperture of 24.0 on the AB system. 

As our primary selection criteria we require: J < 23.7, 
z+ - J > 1.5, J - K s > 0, K s - 4.5^m> 0, no bad 
pixels or cross talk in z + or J bands, and less than a 
2a detection in Bj, g + , and Vj bands. This selection 
removes many galactic stars which have blue J — K s and 
K s — 4.5/im> colors and most low-redshift galaxies 
which are detected in the optical bands. However, our 
sample still contained over 3500 objects, many of which 
appear to be dusty or passive z ~ 2 galaxies or galac- 
tic L and T dwarf starts based on their spectral energy 
distributions (SEDs) (see Section 5). Postage stamps in 
all 22 bands and a weighted average combination of the 
0.3 — 0.8/xm bands were visually inspected to search for 
low level optical flux, resulting in a sample of 22 objects 
with no clear optical detections. Of these 22 objects, 3 
had flat SEDs from 1.2-4.5^m, consistent with z > 7 ob- 
jects. The remaining 19 candidates have very red SEDs 
implying they could be dusty z ~ 2 galaxies. The posi- 
tions of the objects are given in Tableland the F814W, 
and z + limits along with measured fluxes for the 3 objects 
selected for follow-up are tabulated in Table [2] The 2a 
limiting fluxes in a 3" aperture for the other broad band 
images are u* > 27.3, Bj > 27.9, g+ > 27.2, Vj > 27.1, 
r + > 27.3, and i + > 26.9. The intermediate and narrow 



bands have limiting magnitudes ~ 1 — 2 mag brighter 
than the broad bands ( jllbert et al]|2009[ ). 

The color-color plot tor the selected objects are com- 
pared with those from the literatures in Figure [TJ 
Postage stamps of the ground based 0.3-2.2/jm images 
along with Spitzer IRAC Chi (3.6/rni) and Ch2 (4.5/um) 
are shown in Figure ]2] HST-ACS stamps are shown in 
Figure [3| All three objects have Full Width at Half Max 
(FWHM) measurements in the K s band which are larger 
than that measured for nearby stars. However, only ob- 
ject 3 is inconsistent with a star based on the SExtractor 
stellarity index which uses a neural network to consider 
the effects of signal-to-noise and confusion on the size 
measurements. Object 1 is near several foreground, op- 
tically detected sources which partially confuse the IRAC 
photometry for this source. In the 3.6 and 4.5/im bands 
fluxes were recovered using 1.4" diameter apertures cen- 
tered on the K s band position and subsequently correct- 
ing for the missing flux. However, the fluxes at 5.8 and 
8.0/im could not be reliably measured without significant 
de-convolution. Objects 2 and 3 are isolated. Although 
object 2 is technically detected in the ground based i + 
image we included it in our sample because the detection 
is marginal and it is not detected at > 2a in the ACS 
F814W data or the z + band data. 

None of the sources is detected in the XMM 0.5-10 keV 
band X-ray data at a limit of 7.26 x 10 -16 ergs s _1 cm~ 2 . 
Sources 2 and 3 fall into a region with available Chandra 
0.5-10 keV band x-ray imaging and are not detected at 
a 3a upper limit of 7.4 x 10 -16 and 1.0 x 10~ 15 ergs s _1 
cm~ 2 respectively. 

Source 2 is detected in the deep MIPS GO-3 imaging 
of the cosmos field with a flux of 66 ± 12/ xJy at 24/im 



(LeFloc'h et al. 2009 Sanders et al. 20071. Sources 1 



and 3 are undetected in the GU-3 MiPS 24/zm images 
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Object 1 Object 2 Object 3 
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Fig. 6.— The Keck NIRSPEC Nl spectra covering 0.9412- 
1.1097/im is shown with the object positions bracketed by verti- 
cal red lines which are 3.8" apart. The greyscale is set so that 
darker is more positive. For object 1 note the weak line detec- 
tion at 1.0563/^m marked by a red arrow which would place it at 
z = 7.69 if this is Ly-a, in good agreement with the photometric 
redshift. There is no clear detection in the object 2 spectra. For 
object 3 the position of the line seen in the DEIMOS spectra is 
indicated with a blue arrow but there is no detection in the NIR- 
SPEC spectra, however the DEIMOS spectra is ~ 6x as sensitive 
at the NIRSPEC spectra at this position. 

Postage stamps of the Spitzer MIPS 24/im imaging are 
shown in Figure [4] 

Sources 2 and 3 fal l into a region with deep MAMBO 
( Bertoldi et al.|[2007| ) imaging at 1.2mm, and neither is 
detected with fluxes > 2m Jy. Finally, none of the sources 
are detected in the VLA 20 cm data with an RMS of 



16/xJy for object 1 and 10/iJy for objects 2 and 3 ( Schin- 
nerer et al.|2007[ ) or at 90 cm with an RMS of 0.53, 0.42, 
and 0.47 mJy for objects 1, 2, and 3 respectively (Smolcic 
et al. in prep). 

3. SPECTROSCOPIC OBSERVATIONS OF CANDIDATE 
OBJECTS 

Near-IR spectra were obtained on February 22 and 23, 
2008 under photometric conditions with 0. 5" -0.7" seeing 



using NIRSPEC on the Keck-II telescope ( McLean et al 
19981. The data were collected in low-resolution long 
slit mode using the 0.76" slit with the Nl and N5 order 
blocking filters. The resulting spectral resolution was 
R ~ 1500 with a pixel scale ~ 1.6A and ~ 2.7A in the 
Nl and N5 setting respectively and a spacial pixel scale of 
0.18" . The grating was tilted to cover 0.94-1.10/im and 
1.52-1.80/zm in Nl and N5 respectively. These ranges 
were chosen to cover Ly— a (1216A) at 6.7 < z < 8, 
[O II] (3727A) at 1.5 < z < 2, and H-a (6563A) at 
1.3 < z < 1.8. A nearby bright star was used to guide 
the acquisition. For objects 2 and 3 the slit PA was set 
so we could observe the alignment star and the object 




Object 2 Object 3 
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Fig. 7.— The Keck NIRSPEC N5 spectra covering 1.5291- 
1.7953/im is shown with the object positions bracketed by red ver- 
tical lines which are 3.8" apart. The greyscale is set so that darker 
is more positive. For object 1 the object position box is split in half 
with the foreground object on the left, and the target object on the 
right. Note the line in the foreground object spectra at 1.560/im 
indicated by a red arrow. This places the object at z = 1.37 if this 
is H— a. There is no clear detection in the object 1 or 3 spectra. 

simultaneously in at least one of the dither steps. The 
data were taken in an "ABBA" dither pattern with a 
15" throw along the slit. The integration times for ob- 
jects 1, 2, and 3 respectively are 2h, 3h, and 3.75h in Nl 
and lh for all objects in N5. The data were redu ced us- 
ing th e NIR SPEC reduction package described in [Becker] 
et al. (2006). Once reduced, the slits were rectified and 
co-added using a bright object positioned on the slit to 
register the images. 

Latent images from prior observations can be a prob- 
lem in NIRSPEC data. To ensure these were not affect- 
ing the final results we inspected each frame for latent 
images. We also split our data by night, dither pattern, 
and dither position, inspecting each of these stacks for la- 
tent images which appeared in one frame but not others. 
Latent images from a bright standard used to determine 
the slit function were only visible in the first dither Nl 
dither pattern of the spectra taken on object 3 on Feb 22 
and the first Nl exposure taken on object 3 on Feb 23. 
The latents were not at the expected object position. 
Optical spectra were obtained for o bjects 2 and 3 with 



DEIMOS on the Keck-II telescope ( |Faber et ai][2003 l 
on November 23, 2008 under partially cloudy conditions 
and November 25, 2008 under clear conditions with ~ 1" 
seeing. The data were collected with 1.0" wide slits, the 
the 8301/mm grating tilted to 7900A and the OG550 
blocker. The resulting spectral resolution was R ~ 2370 
with a pixel scale of ~ 0.47A a spacial pixel scale of 
0.1185" . The data were reduced with a modified ver- 



sion of the DEEP2 DEIMOS pipeline (Marinoni et al 
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Fig. 8. — The Keck DEIMOS spectra covering 0.64-1.02/xm are shown with the object positions bracketed by red horizontal lines at the 
end of each spectra. The greyscale is set so that darker is more positive. Note the possible line detection in object 3 at 0.9668^m indicated 
by the red arrow and enlarged in a panel at the right of the figure. This suggests a rcdshift of 6.95 if the line is Ly-ct. 



TABLE 3 

Line Sensitivities and un-obscured SFR sensitivities 



ID 


0.6-1/im 




0.95-1.1/jm 


1.5-1.8/Ltm 


SFR-[0 II] (M© yr- 1 ) 


SFR K-a (M Q yr" 1 ) 




erg cm~ 2 s" 


-i 


erg cm" 2 s~ 1 


erg cm" 2 s" 1 


1.0 < z < 2.0, 3.1 < z < 3.8 


1.3 < z < 1.7 


1 






11 - 122 X 10" 18 


1 - 10 X 10~ 17 


1-49, 12-190 


1-16 


2 


1 - 10 X 10" 


-18 


9 - 100 X 10" 18 


1 - 10 X 10" 17 


1-40, 12-190 


1-16 


3 


1 - 10 x 10" 


18 


6 - 70 x 10~ 18 


1 - 10 x 10" 17 


1-28, 12-190 


1-16 



2001 1. In addition to the standard processing this mod- 
ihed pipeline constructs and subtracts a median back- 
ground and accounts for dithering, positional shifts, and 
photometric scaling when combining the spectra. The 
objects were dithered ±3" along the slit to improve the 
sky subtraction and mitigate ghosting (internal reflec- 
tions) from the 8301/mm grating. Bright objects in the 
mask were used to determine the amount of atmospheric 
extinction on November 23, and integrations with less 
than 90% throughput were removed. Each exposure was 
30min for a total of 4.5h of un-extincted integration time. 

The positions of the NIRSPEC and DEIMOS slits are 
shown in Figure [5] and the 2D spectra are shown in Fig- 
ures^ [7J and[Sl We tabulate the range in line sensitivity 
assuming a 10A extraction box in the spectral domain, 
comparable to the expected observed frame line width, 
and the limits this places on the amount of star formation 
in Table E 

We conducted a close inspection of the spectra, includ- 
ing smoothing and binning by up to 32 pixels (~ 115A) 
in the dispersion direction and 10 pixels (1" ) in the slit 
direction. Since the visible spectral features are weak we 
used two methods to estimate their significance. First, 
we created a variance map including read noise, Poisson 
noise, estimated errors in background subtraction, vari- 
ations in depth from bad pixels and dithering, and noise 
correlation due to re-sampling. This map was used to 
estimate the significance of any feature. As a second es- 
timate of significance we used SExtractor to search for 
low signal-to-noise peaks in the 2d spectra with the ob- 
ject and alignment star positions masked. Using this 
catalog of noise peaks, we then calculate the probability 
that a peak of a given flux would fall within ±1" of the 
object position. 

At the expected position for object 1 a weak line is 
visible at 1.0563/im in the Nl spectra. The significance 
of this feature is 2Aa based on the noise map and 2.6c 
based on the noise peak method. If this line is Ly-a it 



places the object at z — 7.69, however the line could also 
be [O II] at z — 1.83. Nothing is seen in the N5 spectra of 
object 1 at the expected position. However, a line is seen 
at 1.560/xm at the position of a nearby optically detected 
foreground object seen in Figure [2] and noted in Section 
2. The noise map yields a significance of 3.8c for this 
feature and no other peaks at this flux level are found in 
the spectra. This would place the neighboring object at 
z = 1.37 if the line is H— a, whic h is the solution fav ored 
by the photometric redshift from Ilbert et al.| ( |2009[ ) , but 
[O II] at z = 3.19 can not be ruled out. [O III] (5007A) 
or H-/3 (486lA) lines at z ~ 2.2 are unlikely for both 
objects because we would expect to see both lines in the 
spectra. However, we note there is significant scatter in 
the observed [O III] and H-/3 line ratios so this redshift 
can not be entirely excluded at our signal-to-noisc. 

Nothing is detectable in the DEIMOS, Nl or N5 spec- 
tra of object 2. 

The DEIMOS spectra of object 3 shows a feature at 
0.9672/im with a 3.4(7 significance based on the noise map 
and no other features with equal significance found in the 
spectra. This suggests a redshift of 6.95 if the line is Ly- 
a, assuming the blue edge is at the line center. However 
this line could also be [O II] at z=1.60 or less likely Ly-/3 
at z=8.43. Nothing is detected in either the Nl or N5 
spectra for object 3. The sensitivity of the Nl spectra is 
~ 6x lower than the DEIMOS spectra at ~ 0.9672/xm so 
the non-detection does not rule out the line seen in the 
DEIMOS spectra. The sensitivity in the N5 spectra is 
also insufficient to detect H— a at z = 1.60 for the typical 
H-q t o [O II] star forming line ratio of 1.77 (Kennicutt 
1998 1 unless the source were obscured by A„ > 7, which 



is significantly larger than the value implied by the broad 
band colors. 

4. SED FITTING ANALYSIS 

To quantify the probability these sources are at z > 7 
and constrain their physical properties, a SED fitting 
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Fig. 9. — The best fit SEDs and redshift probability distributions are shown for all three objects. The best fit galaxy SED is shown in red 
and a secondary galaxy fit, if one exists, shown in orange. The best fit Quasar template is shown in green and the best fit stellar spectra is 
shown in cyan. The joint probability distribution for all SED types is shown in the bottom right corner of each box for regions where the 
probability is greater than 0.01%. Note the lack of any low redshift solution for object 1 and 3 and the peak at z ~ 1.5 for object 2. The 
templates are empirically extended at rest frame wavelengths greater than 2/im. 



TABLE 4 
SED Fitting results 



ID SED Type Redshift e(B-V) 



Mass 



Age 



1 


Emperical 


7 fi+ - 3 
' ' D -0.5 


0.3 






2 


Emperical 


1 - oy -0.02 


0.5 






3 


Emperical 


6.86l°i 


0.3 






1 


BC03 


7 7+1.2 
'•'-0.4 


0.4 


1 x 10 11 


4-30 Myr 


2 


BC03 


1 - OM -0.02 


0.5 


4 x 10 11 


4 Myr 


3 


BC03 


7 n+0.1 
' - z -0.2 


0.1 


2 x 10 11 


450-600 Myr 



analysis w as conducted using the Le Phare codeP°l Both 
empirical dPolletta et al.|2007l|Bu dav ari et al.|2001[ ) and 
theoretical (IrJruzual fc Charlotfl2003[ 2007) templates 
were fit with up to A y = 3 mag ol obscura tion using 
a 



Calzetti et al. ( 1996 ) or Prevot et al. ( 1984b extinction 



law and an additional 2 175 A bump added as a free pa- 



rameter to the extinction law (Ilbert et al. 2009). Stars 
from the |Bixler et ah] ( |1991| ) and |Chabrier et al.j ( |2000 1 
libraries covering O to T stars were also lit to the pho- 
tometry. The optical detection limits were conservatively 
set at 3cr rather than lcr values since our knowledge of 
z > 1 SEDs are limited. The results of the SED fitting 
analysis are shown in Table H] and Figure |9] 

The J band flux is criticalm the SED fitting analysis 
because it constrains the break strength between t he z + 
and J bands, which in turn constrains th e redshift ( Stan 



way et al. 2008b Bouwens et al. 1120081) . However, this 



is mitigated somewhat by the high-quality Spitzer IRAC 
and CFHT H, K s band data which also place strong 
constraints on the SED shape. To test for the effects of 
photometric offsets we introduced a shift of ±0.5 mag- 
nitude in J band and removed it and the IRAC 5.8 and 
8.0/im bands from the fits. 

The best fit for object 1 is an obscured (A„ = 1.2) star 
forming galaxy or un-obscured AGN at 7.1 < z < 8.9, in 

20 http://www.oamp.fr/people/arnouts/LE_PHARE.html 



good agreement with the potential line at z = 7.7. The 
star formation rate based on the SED fit would between 
190-3000 M Q ?/r -1 , depending on the obscuration correc- 
tion, so this source may be detectable in the millimeter 
if it is a star forming galaxy. If instead the source is a 
Type 1 AGN, the x-ray flux limits the black hole mass 
to < 6 x 1O 7 M at z ~ 7 based on an analysis of Type 1 
Quasars in the COSMOS field (Lusso et al., submitted). 

Removing or decrementing the J band magnitude by 
0.5 mag for object 1 produces photo-z solutions between 
0.6 < z < 1.2 with an extinguished star forming galaxies. 
However, all of these fits imply > 3cr detections in the 
i + and F814W^ bands and > 2cr detections in the z + 
band. Furthermore, they suggest star formation rates 
of 30-1000 M^j/r" 1 , which sh ould lead to a MIPS 24/im 



detection (iRieke et al. 120091). As a result these low-z 
solutions are strongly disfavored by the data and it is 
likely that this object is at z > 7. 

For object 2 the low S/N i + and z + band detections 
and 24/xm detection firmly place it at low redshift. The 
best fit SED is a very young, obscured (Ay=1.55) source 
at z = 1.59 with 15 M Q j/r _1 of star formation, in broad 
agreement with the 24/im flux. 

The best fits for object 3 place it at 6.7 < z < 7.3 with 
a 200Myr old relatively un-obscured (A„ = 0.3) stellar 
population or an AGN type spectra. The star formation 
rate implied by the SED fit is ~ 200 M 2/r _1 , comfort- 
ably within the mm detection limits. However the weak 
H band flux relative to the J and if is a quandary. This 
could be explained by an unusually strong 2175A dust 
absorption feature, but, evidence suggests that dust at 
these high redshifts does not have a stron g 2175A feature 
( |Stratta et al.|2007|paiolino et alj2004|) . Alternatively, 
the source could be an AGN at z > 8 with Ly-a and CIV 
falling in the J band and H band covering a region of rel- 
atively line free continuum. Following this interpretation 
the line seen in the DEIMOS spectra would be Ly-/3 at 
z = 8.43, whic h is also seen in deep spectra of the z > 6 
SDSS quasars ( Fan et aT||2006 ), but would imply a very 
strong Ly— a line flux to overcome the strong absorption 
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Fig. 10. — The photometric redshift distribution from Ilbcrt et al. 
1 2009| for z « 7 objects selected by a z+ - J > 1.5, J - K > 0, 
, l s — 4.5/xm> color cut is shown. Note the peak at z ~ 1.8 where 
the 4000A break and dust obscuration is confused with Lyman-a 
forest absorption at z ~ 7. 

from the intergalactic medium at these redshifts. If this 
is an AGN, the Chandra x-ray flux limits the black hole 
mass to < 6x 10 7 M Q at z ~ 8.4 (Lusso et al., submitted). 

Removing the J band data point resulted in a fit at 
2 < z < 4 by an obscured (A„ = 0.3) - 200 Myr old 
stellar population with ~ 200 M ?/r -1 of star formation. 
However, the fits imply > 3er detections in the i + and 
F814W bands and a > 5tr detection in the H band which 
is not seen. 

5. CONTAMINATION FROM Z « 7 

The reliability of photometric redshifts is tightly linked 
to how well the input SEDs represent the galaxy popu- 
lation. Photometric redshifts and photometric selections 
have thus occasionally erroneously identifi ed exotic low- 
redshift objects as high-redshift system s 



et al. 



2i)ii8l |Mobasher et al.|[2005| |Pell6 et al.||2004| |Uickinson 



et al.||2000p . For z > 7 galaxy searches this is especially 
worrisome in view of the large numbers of opticall y faint 
galaxies in the "redshift dese rt" (1 < z < 2.5) ( |Steidel| 
et al.||2004| |Daddi et aT|2004[ ) and late type stars which 
have similar colors and poorly constrained SEDs. These 
contaminating objects are difficult to study for the same 
reasons as z > 7 objects; their strong spectral features 
are shifted to NIR wavelengths where the night sky is 
brighter, the density of night sky lines higher, and instru- 
ments have a smaller fields of view and lower sensitivity. 

To understand the limitations of our selection and the 
SED models we selected a sample of optically bright 
(i + < 26.5) object s from the COSMOS i-band catalog 



(Capak et al.|2007| which meet our z + band dropout se- 

i r-i r-. j — ' . , ~ . . t~ ^ , - ~\ 



5i ■ i ■ 

^ _ Star Forming f ■ 

3- . 



i 2 



z-J> 1.5 

z-I > 1 .5 and Stanway et al. (2008) cuts 
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-1 




Fig. 11. — The (B — z + ), [z — K s ) color-color diagram for an i 
band selected sample is shown. Objects with z + — J > 1.5 colors 
are marked in red, and those detected in all 4 IRAC bands meeting 
the NIR and IRAC color cuts suggested by |Stanway et ajT] ( |2008b] l 
are plotted in green. The position of the sources meeting the z- 
drop criteria indicate they are obscured star forming galaxies and 
passive galaxies at z ~ 2. 



(Figure 12 ) which were formulated to remove passive ob- 
jects at z < 2.5 and late type stars. These cuts succeed 
in removing a large number of objects, but fail to re- 
move a significant fraction of the contaminating popula- 
tion with the steepest spectral slopes. These objects are 
likely to have a steepe r dust extinction curve than the 
locally observed value ( |Siana et al.j2008 ) or be high red- 
shift analogs o f the dusty red ga laxies seen in the A901/2 
super cluster fWolf et al.||2009| ). 

Another possible explanation are objects with extreme 
line equival e nt widths such those presented in Kakazu 
et al.| ( |2007[ ) ; |Hu et al.| ( f2009| . However, at z ~ 1.8 these 
objects would mimick the blue J — H colors tho ught to 
be unique to young stellar popul ations at z > 7 (Oesch 



et al. 2010 Bunker et al. 20101). Furthermore, unlike 



other sources of contamination, these objects are likely to 
be more common at fainter magnitudes because they are 
thought to originate in low-mass, low-metalicity galax- 
ies. However, the same strong emission lines would also 
mean these sources could be easily identified with deep 
spectroscopy, if the sources are bright enough. 

Based on this analysis a range of NIR and IRAC color 
cuts can reduce the number of contaminating objects, 
but not completely remove them. Furthermore, many 
of the z « 7 objects exhibit blue the J — H colors 



The intermediate band data provides a low resolution 


high redshift 


(Oesch et al. 


12010 Bouwens et al. 


2010 


spectra (R~ 20) for these bright objects, allowing us 


Bunker et al. 


2010). Therefore, blue NIR colors 


and a 



shifts as shown in Figure |l0) 
jects have p hoto metric re< 



The majority of these ob- 
ifts in the 1.4 < z < 2.5 



range (Figure 10 ), so a BzK diagram (Figure 11 ) p rovides 



information on the properties of these galaxies (Daddi 



et al.||2004 |. The BzK diagram indicates a mix of ob- 
scured star forming and passive galaxies at z ~ 2. 

Given this contamination we i nvestigated the J, Kg 
and IRAC color cuts proposed by Stanway et al. ( 2008b I 



elude the strongest selection criteria presented here, and 
in the literature, are the non-detections in the i + and 
z + bands which rule out low-z objects based upon the 
strength and steepness of the detected break. To explore 
the necessary depth in i and z bands required to rule out 
low-redshift objects, in Figure [13] we take the photome- 
try for sources meeting our z-dropout criteria and scale 
them to a J band magnitude of zero . This figure shows 
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Fig. 12. — The NIR and IRAC colors of galaxies detected in i 
band with z + — J > 1.5 colors arc plotted in black along with 
our candidate objects in red. The color selections suggested by 
|Stanway et al.| l |2008b|| are shown as cyan lines. Objects removed 
by tne criteria snown in the top panel are marked in green in the 
bottom panel. Note a significant number of objects meet both color 
criteria but are not removed by the suggested color cuts. These are 
likely o bscu red sta r forming and passive galaxies at 1 < z < 3 (See 
Figures [To] and [TT|. 

a significant number of low-z objects have z + — J > 2 
and i + — J > 3.5 colors. So a survey must be at least 3 
magnitudes deeper in z + and 4 magnitudes deeper in i 
than J to adequately select z > 7 galaxies. 

Based on the limits required in i and z bands to rule 
out low redshift objects, most surveys have reached their 
maximum depth for z > 7 searches unless deeper op- 
tical data is obtained. Assuming 2a limits are suffi- 
cient to rule out low-z objects, the COSMOS i band 
data are sufficiently deep to probe to J = 23.7, the 
depth of our present survey. The Hubble Ultra Deep 



Field (UDF) ( |Beckwith et al.|2006[ ) obtains 2a depths of 
F775W=31.1, and F850LP=30.3, which is sufficient to 



remove contaminating objects brighter than Ff60W of 
27.1, well matched to the present NICMOS data. Finally, 
the Great Observatorie s Origins Deep Survey (GOODS) 
(iGiavalisco et al ||2004| F775W and F850LP data reach 
2a depths of 28.9 and 28.6 respectively, sufficient for the 
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Fig. 13. — The (J-A) colors of z << 7 objects selected by a 
z + — J > 1.5, J — K s > 0, K s — 4.5^tm> color cuts are shown. 
Note that the i band data must be 4 magnitudes deeper to remove 
these objects from the sample. 

depth of the pres ent ground based N IR imaging with a 
limit of J = 24.9 (iBouwens et al. 112007). 



Relying on deep i and z band imaging is problematic 
as we probe deeper due to foreground objects which will 
occupy an increasing fraction of the sky. This effectively 
reducing the sample completeness by removing sources 
that are superimposed on faint foreground sources by 
chance. In our ground based COSMOS images we lose 
> 38% of the un-masked sky area to objects detected in 
the ground based optical bands. The situation improves 
in space due to increased resolution and decreasing size of 
faint galaxies. Assuming the size distribution as a func- 
tion of magnitude given in Jouvel et al. ( 2009 ) , using the 
faintest magnitude bin to extrapolate to F775W=31.1 we 
find that > 10% of the UDF area will be lost to objects 
brighter than F775W=3 1 . 1 . 

Therefore, the z — J, J — K, and F850LP - F110W, 
i^HOTU — F160W based z > 7 search techniques used 
from the ground and space may have reached there prac- 
tical limits. Deeper z ~ J color searches with Ultra- 
Vista, the United Kingdom Infrared Deep Sky Sur- 
vey (UKIDDS), or Wide Field Camera 3 (WFC3) will 
likely be dominated by contamination unless significantly 
deeper i and z band data are obtained, and even then 
they are likely to be incomplete due to confusion with 
foreground galaxies. However, narrower a djacent filter 



sets such as that used in the WFC3-UDF (Oesch et al 
2010 Bunker et al.||2010| ) should decrease the amount of 



contamination by better differentiating the smooth roll- 
over in the SED of obscured and old stellar populations 
at z ~ 1.8 and the sharp breaks caused by the Ly— a for- 
est at z > 7. However, objects with strong line emission 
will still be problematic. 

6. DISCUSSION 

6.1. Are these objects at z > 7? 

Based on the present photometric and spectroscopic 
data it is unclear whether or not the objects presented are 
genuinely at z > 7. All three objects meet the color selec- 
tions suggested by previous investigations dStanway et al 



2008b 



et al. 



Oes ch et al.|[2009| |Labbe et aLl|2006[ |Bouwens 



2008 r However, the 24/zm and marginal i + de 
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tection of object 2 strongly suggests it is an obscured 
source at z ~ 1.6. The SEDs of objects 1 and 3 are in- 
consistent with known z < 7 galaxies, AGN or stars in 
our galaxy, and the features seen in the NIRSPEC and 
DEIMOS spectra strongly suggest objects 1 and 3 are at 
z — 7.7 and z — 6.95 respectively. 

The broad band colors and SEDs of these objects are 
indistinguishable from z > 7 candidates found in the 



Bouwens et al. 


2008). 


to those ot Uesch et al. 



2009 



mag 



we estimate a space density of 2.51 ± 1.78 x 10 
Mpc~ 3 for objects with M uv = -23.85 (see Figure 1 141) . 
This space density is consistent with the fain t end counts 
of z ~ 7 galaxies from the WFC3-UDF data (lOesch et al. 



2010 Bunker et al.||20i0 | if the galaxy luminosity iunc- 

the halo mass function shown 
This is significantly higher 



tion follows the shape of 
as a black line in Figure 14 



than the space density implied by fits to smaller area 
data, but not above the theoretical limit for the density 
of these sources. 

If these objects are at z << 7 they point to a popula- 
tion of stars and/or 1.4 < z < 2 galaxies with unexpected 
colors that are contaminating our sample and the objects 
discovered in fainter samples may be contaminants. This 
population of unusual objects must be understood before 
a reliable z > 7 selection can be made. 

We also explored the possibility that object 1 is 
strongly lensed rather than intrinsically bright. To do 
this we constructed a 3D lens model with "Lens Tool" 



( Julio et al. 2007 Kneib et al. 1993 1 using the stellar mass 
and photometric redshhts ot the three nearby foreground 
objects . We assumed a dark-matter to stellar-matter ra- 



tio of 20-1 a nd a Navarro- Frenk- White (NFW) (Navarro 
et al. 19971 profile scaled to the observed light profiles 



of the galaxies. This produced a magnification of ~ 1.2, 
and up to ~ 2 if the galaxies were several times more 
massive than estimated. As a result, we believe lensing 
magnification is an insignificant contribution. 

The lack of conclusive spectra is troubling because 
these objects are at the practical spectroscopic limit for 
8-10m ground based telescopes. Integration times of 
~ 25h (~ 3 nights) per object would be required to reach 
a S/N of 5 for the features seen in these spectra. Multi- 
band NIR ground based spectrographs will allow better 
probes by simultaneously recording spectra from a broad 
wavelength range, however gaps will sill remain in the 
spectral coverage due to atmospheric absorption. Multi- 
object spectrographs will not provide a significant ad- 
vance because the sky surface density of these sources is 
low, even at fainter magntitudes. However, multi-object 
spectrographs will help quantify the contaminating pop- 
ulation and increases in throughput provided by mod- 
ern IR detectors should help significantly. Ultimately, 
ground based 30m telescopes or space based spectroscopy 
from HST or JWST will be required to make significant 
progress towards spectroscopcially confirming this pop- 
ulation of objects. 

Alternatively, a very promising means of confirming 
the redshift is through the [CII] 158/im line. This line is 
the dominant gas cooling line in the interstellar medium, 
tracing the cold neutral medium and photon-dominated 
regions at the interface with active star forming regions, 
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Fig. 14. — Constraints on the z ~ 7 galaxy luminosity function 
from ground and space based data JOesch et al,| 2010 Ouc hi et al.| 
|2009||Oesch et al |2009] |Bouwens et al.liMUS) are shown along with 
the expected di s tribut ion ot dark matter haloes (black line) from 
|Bouwens et al.| (|2008| ) re-normalized to the WFC3 points at the 
taint end ot the luminosity function. The dotted lines indicate the 
limiting magnitude and survey volume limit for the various surveys. 
The error bars on upper limits indicate the point at which 3 objects 
would be detected in the surveyed area to give an estimate of cosmic 
variance. Note the GOODS and UDF data are insufficient to probe 
the bright end of the luminosity function and the density of sources 
is consistent with the expected dark matter halo distribution. 

and it is now being routi nely detected at z > 4 using ex- 
isting sub-mm telescopes (|Maiolino et al.|2009l Ion o et al. 
[2006| |Maiolino et al.||2005| . | Walter et al.| ( [2009| ) show 
that the [CHJ line is a critical tracer of star formation 
in the first galaxies. Bootstrapping from the empirical 
relationship between [CII] luminosity and star formation 
rate based on existing high-z samples, we expect a peak 
line flux density of ~ 1 to 3 mJy for these sources. Such 
a line could be detected with the Plateau de Bure inter- 



ferometer at 240 GHz in about 10 hours (Maiolino et al. 
20091. Furthermore, once operational hi 2013, ALMA 
will allow for high resolution imaging of the gas on sub- 
kpc scales in even short integrations. 

6.2. Implications for galaxy studies of the z > 7 
universe 

A change in the shape of the luminosity function at 
z > 6 is not unexpected. In Cold Dark Matter (CDM) 
models the knee of the galaxy luminosity function is usu- 
ally explained by feedback which truncates star forma- 
tion in the highest mass haloes. If feedback is inefficient 
in the first generation of galaxies, their luminosity func- 
tion should flatten out to follow the power-law halo mass 
function at these redshifts. The galaxy number density 
meas ured in this paper is below the expected ha l o den- 
sity, (fC ooray & Ouchi 2006; Bo uwens et al.||2008j |Ouchi 
et al.|2TjL)9l |(Jesch et al.||2010 ), so we can not rule these 
objects out based on space density alone. The most likely 
feedback mechanisms are Active Galactic Nuclei (AGN), 
supernova, and stellar winds. But, in the early universe 
these mechanisms are likely to be inefficient. AGN feed- 
back does not become efficient until the central black 
hole grows to a large enough mass to provide significant 
energy output. Simulation suggest this takes ~200-500 
Myr, lea ving plenty of time for bright galaxies to exist 
at z > 6 ( Li et al.| 2007 ) . Furthermore, simulations have 
shown that supernova do not contribute enough energy 
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to truncate star formation, and may act ually enhance 



star f ormation in a gas rich environment (Wise & Abel 
2008 1. Finally, gas in the early universe will have very 



low metalicity, making winds less efficient at removing 
material. 

Several empirical studies also support the existence of 
bright starbursts at z > 7. Studies of the z ^ 4 — 6 mass 
function point to a epoch of star for mation at z > 6 
( [Stark et al.1|2009| |McLure et aL||2009| ), and fe w studies 



suggest even more massive ob jects may exist (Wiklind 



et al.||2008[ Eyles et al. 2007), but there is significant 
disagreement on th is latter point (see Dunlop et al. 2007 
Chary et al. 2007). Further more, z ~ 6 qua sars also 



exhibit signific ant enrichment ( Ku rk et al. 2007 ) and high 
dust contents ( Riechers et al.|20"07[ ) implying a significant 
earlier generation of stars was present. 
If we take the density of M~ 10 n M^ galaxies to be 



1 x 10~ 6 dex^M joc- 3 at z ~ 6 (Stark et al. 



2009 



McLure et al. 2009), the continuos star formation rate 
required to form these objects between z ~ 20 and z ~ 6 
is > 130 M Q yr _1 , which implies an intrinsic J band 
magnitude of J < 23.6 at z = 7, placing an upper limit on 
the density of these systems of ~ 4 x 10~ 7 mag -1 Mpc -3 . 
Star formation in bursts would make these objects even 
more luminous, but also rarer because they are short 
lived. Assuming a minimum formation time of 100 Myr 
for a 10 11 M Q galaxy yields a lower limit on the space 
density of ~ 5 x 10~ 8 mag" 1 Mpc -3 . Our measured 
value is reasonably between these two extremes. 

In addition, there is significant disagreement about the 
evolution of the bright end of the luminosity function at 
z > 4, so the scale of disagreement between our result 
and those at z ~ 6 is unclear. Deep ACS surveys sug- 
gest the density of br ight galaxies is dropping at z > 4 
(Bouwens et al. 2008) . However, these surveys cover a 
very small volume and hence lack significant constraints 
on the bright end of the luminosity function. Larger 
area ground based surveys prov ide mixed conclusions 
with |Sawicki fe Thompson] ( |2006[ ) and jlwata etal] ( |2007 1 
point ing to no evolution m the density of brigh t galaxies 
while |Yoshida_^aI1 fl2006[ ) and |Capak| (|2008| ) find mild 
evolution. McLure 'et al.| ( |2009p find strong evolution, 
however their results are based on photometric redshfits 
which can miss a significant number of bright objects at 
these redshifts (Capak et al. in prep). In any case, the 
continuum luminosity function of spectroscopically con- 
firmed Ly-a emitters shown in Figure [15] also places a 
lower limit on the luminosity function at z ~ 5.7 which 
is higher than the estimate from many studies at z ~ 6 
as is also noted by Hu & Cowie (2006). 

While it is numerically possible to fit a luminosity func- 
tion to our data at z ~ 7 the results would be mislead- 
ing given the present uncertainty in the data. In ad- 
dition to the uncertainty in the redshift of our sources, 
other surveys are likely dominated by cosmic s hot noise 



at the brigh t end due to their small volumes (Overzier 
et al. 2009). At z > 5 the luminous galaxy popula- 
tion is extremely clustered, with bias facto rs of ~ 6 — 10 
( |Hildebrandt et al.|2009 



Ouchi et al. 12004). This means 



the mode (most common) luminosity function probed by 
small area surveys is likely to systematically underesti- 
mate the density of luminosity galaxies because the mas- 
sive haloes which host these massive (bright) systems are 
unlikely to occur in the survey volume. To account for 
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Fig. 15. — Data on the z ~ 6 UV (1500A) luminosity function 
from the literature is shown with LBG samples plotted in black 
| |McLure et al.||2009| [Bouwens et al.|20 07 Shimasaku ct al j 2005 
and z ~ 5.Y LAID samples plotted in red | |rlu (Jowic 2U Uo| [shr 
|masaku et al.||2006[ |Ouchi et al.|[2~008| |. The LAD samples are 
truncated at the continuum magnitude where they begin to be in- 
comp lete. The best fit lum inosity functions at z ~ 4 and z ~ 6 
from iBouwens et al.l (|2007| are shown. We note that the LBG 
and LAID data arc statistically consistent with the exception of 
the |McLure" ct al. ( 2009| results which are based on photometric 
rcdshitts. 

this uncertainty we have set the upper limits in Figure 



14 to the point where they would detect 3 galaxies in 
an attempt to account for the presently unknown galaxy 
clustering properties. 

7. CONCLUSIONS 

The key to any successful Lyman break selection cri- 
teria is differentiating between the Lyman break and the 
4000A break at lower redshifts and separating these ob- 
jects from stars. For a 2 < z < 4 LBG selection the low 
redshift contamination is from z < 0.5 where the galaxy 
population is well understood along with O, B, and A 
stars which are relatively rare ( |Steidel et al.|200"3 ) . For a 
4 < z < 6 selection contamination is from 0.5 < z < 1.5 
galaxies where the population is not as well un derstood 
and M stars which are much m ore common (Stanway 

Furthermore, spectro- 



et al. 


2008a 


Capak et al. 


2004) 



sources found by a 4 < z < 6 LBG selection is more 
difficult. Since these selection effects are not well under- 
stood this becomes a likely source of the large discrep- 
ancies seen between z ~ 5 — 6 galaxy studies (|Stanway| 

" 'Iwata et aTT : 2007). 



et al. 


2008a 


Bouwens et al. 


2007 



1.5 < z < 2.5 galaxies and L,T, and Y dwarfs, all of 
which are poorly understood and extremely difficult to 
spectroscopically confirm. This is mitigated somewhat 
by Spitzer IRAC and K s band data which place strong 



constraints on th e SED shape (ptanway et al. 2008b 



Labbe et al. 2006). However, the strength of these con 



stramts is limited by the lack of a representative set of 
object templates at 1.5 < z < 2.5 and the large uncer- 
tainties on deep Spitzer photometry due to confusion and 
under-sampling. 
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Of our three candidates, we find object 1 is likely, but 
not conclusively at z = 7.7, object 2 is at z <~ 1.6 
and object 3 is very likely at 2 = 6.95. However, the 
faint H band flux of object 3 is puzzling, suggesting a 
strong 2175A dust absorption feature at z = 6.95 or a 
Quasar type spectra with strong Ly-a and CIV emission 
at z = 8.43. The spectral features observed in these 
sources are weak but known galaxy SEDs fail to produce 
a z « 7 solution for these objects, and no objects at 
brighter magnitudes exhibit a similar SED, so together 
the evidence is reasonably strong. 

Whether or not these sources are at z > 7, the dif- 
ficulty in confirming even these bright objects suggests 
that deep small area surveys with small numbers of fil- 
ters provide poor constraints on the z > 7 galaxy pop- 
ulation. If the presented galaxies are at z > 7, then 
galaxies likely follow the dark matter halo mass function 
without a cut-off at the high-mass end as is observed at 
low redshift. As a result pencil beam surveys will fail to 
provide an accurate census of the galaxy population due 
to a lack of volume which prevents them from finding 
the most luminous galaxies and cosmic variance which 
increases their error. If we assume these galaxies are 
contaminated by z << 7 objects, the z-dropout selection 
function needs to be quantified by spectroscopy before 
color selections can be used. If high equivalent width 
line emitters are the main source of contamination, then 
current NIR spectroscopy instruments may be able to 
determine this. However, direct spectroscopic confirma- 
tion of the z > 7 galaxy population will not be possible 
until thirty meter telescopes, JWST, or possibly ALMA 
is available. Either way, medium deep, pan-chromatic 
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